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Abstract 


This  report  presents  the  results  of  research  project  with  the  aim  of  evaluating  the  efficiency 
of  energy  deposition  by  combustion  upstream  blunt  bodies  in  supersonic  Mach  number  2.0-5. 0 
air  streams.  At  these  conditions,  considerable  wave  drag  reduction  with  a  high  efficiency  of  the 
energy  utilization  has  been  demonstrated  recently  experimentally  (see  [2]  and  in  the  modeling 
(see  [2],  [3],  [4]).  In  the  most  discussions,  the  energy  source  is  supposed  to  be  attributed  to 
the  focused  optical  (laser)  or  microwave  discharges  with  an  important  role  of  charged  particles 
produced  by  discharges.  For  example,  it  has  been  claimed  in  [lljthatthepresenceof  only  one  (!?) 
electron  per  10^  neutral  particles  can  lead  to  substantial  reduction  of  both  viscous  and  pressure 
drag. 

The  objective  of  this  study  is  to  assess  the  feasibility  of  utilizing  the  thermal  energy  of 
counterflow  combustion  realized  using  a  blunt  body  with  a  smaller  diameter  pipe  working  as 
the  aerodynamic  spike  injecting  combusti bl e  gases  upstream  the  supersonic  flow.  At  these  con¬ 
ditions,  the  charged  particles  are  not  formed,  and,  drag  reduction,  if  confirmed,  could  be  proven 
not  to  be  the  plasma  specific  phenomenon. 

First  of  all,  the  experimental  evidences  of  drag  reduction  by  counterflow  combustion  in  su¬ 
personic  stream  are  presented.  The  procedure  and  hardware  used  in  experiments  is  described  as 
well.  Then,  to  perform  the  problem  numerical  study  based  on  the  solution  of  full  Navier-Stokes 
equations,  the  Kl  VA-3  reactive  flow  computer  code  has  been  modified  to  simulate  the  counter¬ 
flow  combustion  in  the  M  =2.0  supersonic  stream  using  the  detailed  chemistry  approach.  In  this 
study,  both  experimental  and  theoretical,  the  hydrogen/air  mixture  has  been  analyzed.  The  re¬ 
sults  illustrate  the  considerable  effect  of  the  counterflow  combustion  on  a  flow  structure  before 
the  blunt  body  in  the  M  =2.0  supersonic  stream.  Flowever,  the  intensity  of  heat  release  appeared 
not  to  be  sufficient  enough  to  generate  a  more  effective  "air  spike"  flow  structure  typical  of  appli¬ 
cation  of  the  powerful  optical  (laser)  discharge.  To  increase  the  energy  release  rate,  the  kinetics 
of  the  gasified  boron-based  hydrocarbon  fuel  composition  has  been  analyzed.  The  chemical  me¬ 
chanism  has  been  constructed,  and  its  more  high  combustion  intensity  is  demonstrated.  The 
important  role  of  a  conical  tip  of  the  spike  for  flame  stabilization  has  been  also  demonstrated. 

The  direct  measurement  data  of  drag  reduction  by  counterflow  hydrogen  combustion  in  M=2.0- 
2.5  air  streams  are  presented  and  discussed  together  with  the  numerical  predictions.  These  re¬ 
sults  confirm  that  the  "flame  spike"  could  bean  effective  means  of  drag  reduction  of  blunt  bodies 
in  supersonic  flows. 


1.1  Introduction 


The  concept  of  structural  aerodynamic  spike  is  an  effective  means  to  reduce  the  pressure 
drag  of  blunt  bodies  in  supersonic  stream.  The  example  of  the  concept  realization  is,  e.g.,  the 
Trident  D-5  missilewhich  flying  distance  was  thus  increased  by  550  km.  Other  applications  are 
also  conceivable. 

I  n  the  presence  of  the  structural  spike,  the  flow  structure  is  characterized  by  the  presence 
of  the  oblique  shock  and  the  recirculation  zone  forming  an  aerodynamic  body  profile  similar  to 
that  of  a  contoured  body.  If  the  spike  has  a  non-optimal  length  (e.g.  too  long),  the  truncated 
recirculation  zone  forms  (see  Fig.1.1),  and  as  far  as  the  oblique  shock  is  stabilized  above  the 
front  boundary  of  this  zone,  drag  reduction  effect  will  be  moderated. 


(a)  (b) 


Figure  1.1:  Supersonic  flow.  Moo  =2.0,  structures  beforeablunt  body  with  a  structural  spike.  Two 
types  (depending  on  the  spike  length)  of  flow  separation  has  been  observed:  a)  the  leading-edge 
and  b)  delayed  separations. 


Drag  reduction  is  more  related  to  the  problem  of  construction  of  a  body  surface  of  minimum 
drag.  Among  possible  analytical  solutions,  in  aclassof  axysimmetricbodies,  there  is  the  solution 
with  a  singularity  of  the  body  surface  that  is,  in  fact,  a  structural  spike.  Due  to  a  formal  simi¬ 
larity  between  momentum,  mass,  and  energy  conservations,  it  is  possibleto  show  that  the  drag 
reduction  effect  can  be  achieved  by  varying  momentum,  mass,  and  energy,  i.e.  mass  or  energy 
"spikes"  are  conceivable.  The  effect  is  dependent  on  the  type  and  size  of  non-uniformity  of  the 
flow  around  the  body.  The  wake-like  non-uniformity,  produced,  for  example,  by  the  structural 
spike  due  to  non-slip  velocity  conditions,  forms  stable  flows  with  the  front  recirculation  zones 
that  reduces  drag.  On  the  contrary,  the  jet-type  non-uniformity  (produced,  if  the  body  has  a 
recess  or  a  blind  hoi  eon  the  front  surface)  leads  to  non-steady  oscillating  flow  regimes  used  in 
acoustic  generators/  resonant  tubes. 

Some  examples  of  such  a  generalized  spike  concept  applications  are  presented  below.  In 
Fig.l.2a,  the  effect  of  the  structural  spike  is  simulated  by  a  bead  projected  upstream  the  su¬ 
personic  M=2.0  flow.  The  sub-sonic  wake  behind  the  bid  effectively  forms  the  wake-like  non¬ 
uniformity  of  the  velocity  profile.  The  effect  of  such  a  "spike"  becomes  well  pronounced  when  the 
bid  reaches  the  distance  similar  to  the  "optimal"  length  of  the  spike.  The  injection  of  mass  can 
result  in  the  same  effect.  I  n  Fig.l.2b,  the  drag  reduction  is  illustrated  by  Shkval,  a  cone-shaped 
rocket-torpedo  that  travels  within  a  supercavitating  bubble  flow  produced  with  the  help  of  its 
own  engine  exhaust. 

Thefurther  development  of  the  idea  haslead  recently  to  the  "air-spike"  concept  implemented 
in  two  different  ways  by  Myrabo  [1]  and  Tretjakov,  et  al.  [2].  The  air  spike  can  be  formed  by 
concentrated  energy  (an  electric  arc  plasma  torch  in  the  Myrabo's  case,  a  repetitive-pulse  laser 
beam  in  the  case  of  Tretjakov  et  al.,  see  Fig.l.2c)  projected  forward  off  a  moving  body  producing 
a  "tunnel"  of  low  density,  reduced  pressure  hot  air  in  the  shape  of  a  paraboloid  of  revolution. 
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(a) 


(b)  (c) 

Figure  1.2:  a)  Simulation  of  the  structural  spike  effect  by  projecting  a  bid  upstream  a  superso¬ 
nic  M=2.0  flow,  b)  Drag  is  reduced  by  creating  a  local  "envelope"  of  supercavitating  bubbles  of 
gaseous  combustion  products  in  which  the  weapon  exceeds  a  speed  230  mph  -  a  "bubble"  spike, 
c)  Schlieren  image  illustrating  Tret] akov  et  al.  air-spike  concept  realization  by  concentrating  a 
laser  energy  before  the  body  i  n  the  M  =2.0  flow. 


Such  a  "spike"  has  important  advantages  over  a  structural  spike  due  to  the  fact  that  air  density 
behind  the  blast  wave  is  lower  than  that  behind  the  shock  wave.  The  level  of  the  drag  reduction 
in  the  supersonic  M  =2.0  flow  can  reach  ~  50%  of  the  baseline  drag.  The  effect  depends  also  on 
the  relative  "hot  spot"  position  with  respect  to  the  body,  the  "hot"  spot  size,  and  M  number  of  the 
free  stream. 

Counterflow  combustion,  as  proven  in  [5]  and  illustrated  in  Fig.l.3a-c,  tends  to  weaken  (or 
nearly  suppress)  the  shocks  near  the  blunt  body  face,  thus  reducing  the  drag.  In  this  case, 
the  spike  is  used  to  be  an  injector  through  which  the  combust! bl e  gas  is  injected  upstream  the 
supersonic  flow.  The  positive  effect  of  such  a  "flame  spike"  has  been  validated  both  theoretically 
and  experimentally  and  some  measurement  data  are  presented  in  Fig.l.4for  different  (including 
projectile)  geometries. 

The  combined  effect  of  the  structural  and  "flame"  spikes  is  well  illustrated  in  Fig.l.4a,  in 
which  one  can  see  that  the  drag  reduction  in  the  presence  of  counterflow  combustion  is  better 


(a)  (b)  (c) 

Figure  1.3:  The  M=2.0  flow  structures:  a)  in  the  absence  of  counterflow  injection,  b)  in  the 
presence  of  inert  gas  injection,  (5/^2=0.002,  and  c)  in  the  presence  of  counterflow  combustion, 
(5 i/2  =0.002.  The  hydrogen  mass  flow  rate  amounts  to  0.2%  of  the  air  mass  flow  rate  through  the 
body  cross  section. 
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Figure  1.4:  a)  The  pressure  drag  vs  the  structural  spike  length  for  the  blunt  cylinder  at  M=2.0 
flow:  1-  pressure  drag,  2-  base  drag,  3-  total  drag  with  hydrogen  combustion.  The  pressure  drag 
as  a  function  of  a  relative  mass  flow  rate,  G,  for  fuels  with  different  heat  release  in  combustion, 
b)  for  the  model  geometry:  1-  iy„=120  MJ  /kg,  2-  iy„=16.7  MJ  /kg,  3-  iy„=13.4  MJ ,  c)  for  real 
rounds:  1-  iy„=16.7  MJ  /kg,  2-  i7„=13.6  MJ  /kg,  3-  i7„=9.0  MJ .  Low  heat  values  correspond  to  solid 
propellants  used  in  the  nose  gas  generator. 


pronounced  and  I  ess  susceptible  to  change  in  the  case,  when  the  structural  spike  length  becomes 
non-optimal.  TheFig.l.4b  illustrates  that  the  drag  reduction  effect  can  be  achieved  using  very 
low  mass  flow  rates  of  combustible  gases,  and  is  better  pronounced  for  higher  energy  density 
materials.  In  all  illustrations,  the  relative  mass  flow  rate  is  the  ratio  of  the  injected  fuel  mass 
flow  rate  to  the  air  mass  flow  rate  through  the  body  cross  sect!  on  area.  As  fol  I  ows  from  Fig.l.4a, 
the  concept  proposed  seems  to  be  more  effici  ent  than  the  wel  I  -known  drag  reducti  on  means  usi  ng 
base  burning.  A  wide  range  of  drag  variation  opens  perspective  to  realize  a  guided  flight  using 
low  fuel  mass  flow  rate  counterflow  combustion. 

The  flame  spike,  if  applied  to  a  supersonic  cruise  missile,  could  reduce  drag,  thus,  preventing 
deceleration  of  the  vehicle  driven  by  a  relatively  low  power  engine  with  an  additional  benefit  of 
the  spike  "products"  can  envelope  the  missile  in  a  film  of  slightly  ionized  gas  which  would  be 
impervious  to  radar  pulses,  thereby  rendering  it  electronically  "invisible".  The  infra-red  penalty 
at  a  close  range  would  be  significantly  offset  by  the  high  speed  of  the  flight.  The  missile  could  be 
initially  accelerated  using  a  conventional  solid  propellant  booster. 

1.2  Description  of  Work 

The  purpose  of  the  project  was  to  achi  eve  as  fol  I  ows: 

•  Obtain  via  direct  measurements  the  data  confirming  the  ability  of  counterflow  combustion 
to  reduce  drag  of  blunt  bodies  in  supersonic  M  =2. 0-2.5  air  streams, 

•  Present  the  theoretical  background  to  evaluate  the  drag  reduction  effect  of  counterflow 
combustion, 

•  Determination  (experimentally  or  theoretically)  of  stableenergy  release  conditions  in  high 
speed  flows, 

•  Present  information  available  and  search  for  new,  more  effective  than  hydrogen,  fuel  com¬ 
positions, 

•  Study  of  combustion  completeness  and  possible  positive  effect  optimization. 
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1.3  Experimental  Study 

The  M  =2.0  wind  tunnel  of  thel  nstitutefor  Pure  and  Applied  Mechanics  of  Siberian  Branch  of 
Russian  Academy  of  Sciences,  Novosibirsk  has  been  used  in  the  flame  spike  concept  validation. 
The  detailed  description  of  the  hardware  can  be  found  elsewhere,  and  a  brief  description  of  the 
particular  experiments  is  given  below. 

In  Fig.  1.5a,  it  is  shown  the  experiment  scheme  enabling  one  to  study  the  flow  structure 
in  the  presence  of  the  structural  spike,  the  hysteresis  phenomena-  the  transition  from  one  flow 
separation  pattern  to  another  one  when  the  spike  length  is  varied,  as  well  as  the  effects  of 
the  mass  supply  through  the  spike  and  combustion  on  the  flow  structure  and  drag.  The  spike 
hysteresis  effect  was  specially  investigate  ddue  to  a  possible  consequence  for  the  moving  body 
with  thespiketo  have  different  drag  depending  on  whether  the  body  accelerates  or  decelerates. 


(b) 

Figure  1.5:  a)  Model  in  a  test  section  of  the  wind  tunnel.  Flow  conditions:  M  =2. 0-2.5,  Re«lO^. 
Notations:  1-  model,  2-  aerodynamic  spike,  3-  model  aft  mounting  sting,  4-  pilot  ignition  device, 
5-  optical  window,  6-  reacting  gas  supply,  7-  towards  the  mechanism  of  a  vertical  displacement 
of  the  pilot  ignitor,  8-  towards  the  mechanism  of  a  horizontal  displacement  of  the  spike,  b) 
Aerodynamic  spikes  of  different  shapes:!-  sharp  end  spike,  2-  pipe,  3-  pipe  with  a  conictip. 


The  aerodynami  c  spi  ke  can  be  moved  al  ong  the  model  axi  s  wi  th  the  hel  p  a  coordi  nate  appi  i  ance- 


4 


the  mechanism  for  a  horizontal  displacement.  The  tested  aerodynamic  spikes  of  different  shape 
are  presented  in  Fig.l.Sb.  The  spike  supplying  gases  may  be  argon,  hydrogen  and  air.  The  ex¬ 
periments  with  hydrogen  injection  followed  by  combustion  were  carried  out  as  well.  Pressure 
was  measured  on  the  face  surface  of  the  model  with  the  help  of  inductive  gauges  in  locations 
along  the  radius  and  recorded  by  the  potentiometer.  The  spike  gas  flow  rate  was  measured  by 
thestandard  gauge  based  on  the  diaphragm  pressure  drop.  The  displacement  speed  of  the  spike 
accounted  for  approximately  1  mm/s.  The  flow  regime  change  was  registered  with  the  help  of 
optical  observation  of  the  flow  structures  through  the  windows  using  the  Tdpler's  device  as  well 
as  by  pressure  measurements  along  the  model  face  surface. 

To  initiate  counterflow  combustion,  a  microramjet  pilot  device  (seeFig.l.6a)  has  been  used  as 
a  source  of  high  temperature  products  of  hydrogen-air  combustion.  This  device  operated  reliably 
within  the  range  of  M  =2. 0-2.5.  The  rich  flame  quenching  regime  did  not  take  place  within  the 
limits  of  maximal  possi  bl  e  hydrogen  mass  flow  rate  through  the  igniter  (~  0.2  g/s)  determined  by 
the  design  peculiarities.  Fig.l.Ob  illustrates  the  ignition  device  operation.  Ignition  of  the  main 
hydrogen  injected  from  the  spike  was  occured  by  means  of  the  microramjet  device  movement 
towards  the  hydrogen  injection  position.  After  hydrogen  ignition,  the  igniter  was  removed. 


Figure  1.6:  a)  Schematic  of  the  pilot  ignition  device  (1-  insulator,  2-  ignition  spark  electrodes) 
and  b)  visualization  of  its  operation. 


A  few  technical  realizations  for  combustion  initiation  were  considered.  They  were  based  on 
the  application  of  small  pyrotechnic  ignition  cartridge  (using  silane,  SiH4)  and  a  special  realiza¬ 
tion  (see  [5])  of  the  resonance-igniter  concept  for  gas  heating  in  a  closed  end  duct  similar  to  the 
process  in  the  resonant  tubes. 

Flydrogen  wasthefuel  of  choice  because  of  its  high  energy  release  rate  and  because  of  hydro¬ 
gen  can  be  formed  in  a  considerable  amount  in  the  pyrolysis  of  practical  hydrocarbon  fuels  (see, 
e.g.,  [7]).  First  of  all,  the  positive  effect  of  counterflow  hydrogen  combustion  on  drag  reduction  in 
supersonic  M=2.0-2.5  flows  has  been  directly  measured  and  the  data  aresummarized  in  Fig.1.7. 
F rom  these  data,  one  can  conci  ude  that  drag  of  the  bl  unt  body  with  the  structural  spi  ke  can  be 
reduced  to  60  %  of  the  baseline  drag  when  using  injection  and  combustion  of  asmall  amount  of 
hydrogen.  The  relative  drag  reduction  in  Fig.1.7  is  presented  as  a  function  of  a  parameter  which 
is  a  product  of  the  relative  mass  flow  rateG  and  the  relative  heat  addition  Ho/CpTo,  where  i7o 
is  the  combustion  heat  release.  The  positive  effect  of  drag  reduction  by  counterflow  combustion 
is  comparable  with  the  effect  achieved  with  the  help  of  optical  (laser)  discharge  (see  [8])  for  the 
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Figure  1.7:  Drag  reduction  effect  of  counterflow  hydrogen  combustion  in  supersonic  M=2.0- 
2.5  streams  for  blunt  bodies  of  different  shapes  (cylinder  and  hemisphere+ cylinder)  with  the 
structurlal  spikes  used  asthefuel  injector. 


same  flow  conditions. 

Then,  a  series  of  experiments  was  carried  out  that  enabled  one  to  determine  the  mini¬ 
mum  sizes  of  the  blunt  body  with  the  aerodynamic  spike  that  makes  it  possible  to  realize  the 
steady-state  hydrogen  combustion  in  air  stream  with  parameters  Moo  =2. 0-2.5,  Ptot  =2.5-45  atm, 
Ttot  =300  K.  The  results  are  summarized  in  Fig.1.8.  Fig.l.Sa  illustrates  the  relationship  between 


(a)  (b) 

Figure  1.8:  Regions  (in  the d-L  and  Pq-G  planes)  of  hydrogen  steady  counterflow  combustion 
in  the  supersonic.  Moo  =2.0  stream:  a)  the  minimum  diameter  of  the  body,  and  the  spike  length; 
b)  the  hydrogen  mass-flow  rate  and  the  total  pressure.  Boundaries  of  the  flame  blow-off:  1  -  a 
"lean"  flame  blow-off,  (A  >1.0),  2-  a  "rich"  flame  blow-off,  (A  <1.0),  3-  a  "lean"  flame  blow-off  for 
a  free  counterflow  jet. 


the  body  diameter,  d,  and  the  dimensionless  length  of  the  structural  spike,  I  =  L/d,  where  L  is 
the  spike  length,  in  the  region  of  stable  combustion.  The  boundary  3  corresponds  to  the  flame 
which  is  stabilized  by  a  separation  zone  formed  by  a  free  hydrogen  jet  not  connected  with  the 
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separation  zone  in  front  of  the  blunt  body  with  the  spike.  The  "lean"  flame  blow-off  takes  place 
when  the  spike  is  longer  than  the  optimal  one.  Fig.  1.8b  shows  the  regions  of  stable  combustion 
in  the  G -p  plane:  I  is  the  boundary  of  the  stable  counterflow  combustion  for  a  free  jet,  1 1  isthe 
stabilization  region  for  the  blunt  body  with  the  spike.  This  region  has  two  boundaries  marked 
by  land  2.  Theflow  structures  corresponding  to  counterflow  hydrogen  injection  and  combustion 
are  presented  in  Figs  1.9-1.10. 


(a)  (b) 

Figure  1.9:  a)  Gasdynamic  M=2.0  flow  structure  for  the  inert  cointerflow  jet:  1-  bow  shock,  2- 
contact  surface,  3-  shock  of  the  jet,  4-  region  of  flow  expansion,  5-  re-attachment  shock,  6-  flow 
separation  zone,  8-  injector,  b)  Dependence  of  the  main  geometrical  parameters  of  the  inert  jet 
on  the  gas  injection  pressure. 


Data  presented  in  Fig.  1.9  correspond  to  the  case  of  the  inert  jet,  in  Fig.  1.10  -  to  the  case  of 
counterflow  combustion.  These  data  were  obtained  with  the  help  of  schlieren  photographs  of  the 
fluid  structure.  Geometrical  characteristics  of  both,  inert  and  reacting  counterflow  jets  are  well 
generalized  using  the  dimensionless  parameter:  a  ratio  of  the  total  pressure  in  the  jet,  Pgj  to 
the  pressure  behind  the  normal  shock  at  the  parameters  of  the  free  stream,  Pg'.  The  effect  of 
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(a)  (b) 

Figure  1.10:  a)  Gasdynamic  M  =2.0  flow  structure  for  the  reacting  cointerflow  jet:  1-  bow  shock, 
2-  contact  surface,  3-  shock  of  the  jet,  6-  flow  separation  zone,  7-  flame,  8-  injector,  b)  Dependence 
of  the  main  geometrical  parameters  of  the  inert  jet  on  the  gas  injection  pressure. 
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combustion  is  well  pronounced:  for  this  case,  the  distance  to  the  normal  shock  3  and  the  contact 
surface  radius  2  are  longer  than  for  the  inert  case  due  to  the  lower  density  in  the  recirculation 
zone  of  the  reactive  flow.  This  also  means  that  the  shock  wave  inclination  for  the  case  of  the 
reactive  flow  becomes  less  that  testifies  to  drag  reduction  caused  by  counterflow  combustion. 

Thus,  one  can  conclude  that  the  possibility  of  drag  reduction  using  the  counterflow  combus¬ 
tion  in  supersonic  streams  realized  in  a  form  of  the  "flame"  spike  is  confirmed  experimentally. 
The  size  of  the  body  and  and  the  spike  length  can  be  selected  to  guarantee  stable  combustion 
for  given  free  stream  conditions.  The  effects  of  the  mass  flow  rate  and  the  energy  deposition 
variations  for  the  model  fuel  (hydrogen)  are  studied.  It  has  been  demonstrated  that  drag  of  the 
blunt  body  with  the  "flame"  spike  in  M  =2.0-2.5  streams  can  be  reduced  up  to  60  %  of  drag  of  the 
baseline  blunt  body  with  the  structural  spike.  At  that,  the  blunt  body  shape  does  not  play  an 


(a) 


(b) 


(d)  (e) 

Figure  1.11:  Spike  hysteresis-  M  =2.0  flow  structures  and  pressure  diagram.  Flow  structures  for 
a)  the  optimal  length  structural  spike,  b)  the  long  structural  spike,  d)  the  short  length  flame 
spike  (G  =0.015),  e)  the  long  flame  spike.  In  c)  it  is  presented  the  total  pressure  vs  spike  length 
(in  mm)  record. 
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essential  role.  The  interesting  feature  of  the  experimental  data  obtained  is  the  effect  modera¬ 
tion  (see  Fig.1.7),  if  the  fuel  amount  injected  or  the  amount  of  energy  released  will  exceed  some 
specific  values.  Thiseffect  will  be  discussed  later  on. 

The  spi  ke  hysteresi  s  effect  I  s  very  i  mportant  for  practi  cal  appi  i  cati  on  due  to  the  consequence 
that  the  moving  body  with  the  spike  may  have  different  drag  depending  on  whether  the  body 
accelerates  or  decelerates.  I  nformation  presented  in  Fig.  1.11  can  be  used  to  clarify  the  matter. 
Hysteresis  was  shown  to  exist  between  the  two  types  of  flow  separation  caused  by  the  presence 
of  the  spike:  the  first  one  occurs  (see  Fig.  1.11a)  when  the  spike  has  the  optimal  length,  the 
second  one  ((see  Fig.  1.11b)  is  typical  of  the  long  spike.  Due  to  the  hysteresis  (see  Fig.  1.11c), 
it  is  possible  to  enter  the  high  drag  state  from  low  drag  at  the  same  spike  length,  although 
the  reverse  is  not  true.  Date  presented  in  Fig.  l.lld-e  are  characteristic  of  the  counterflow 
combustion  reducing  this  effect  substantially:  for  two  flame  spikes  of  different  lengths  this  is 
no  transition  from  one  flow  separation  regime  to  another  observe.  This  explains  the  results 
presented  in  Fig.  1.4and  illustrates  the  uniquefeature  of  the  flame  spike.. 


1.4  Theoretical  Study 

At  first,  the  counterflow  combustion  in  a  supersonic  flow  was  studied  numerically  with  the 
hel  p  of  a  computer  code  descri bed  in  [2].  I  n  this  case,  the  effect  of  combustion  has  been  simulated 
by  introduction  of  the  heat  source  term  in  the  energy  equation.  The  typical  results  obtained  by 
applyingthecodetoM  =2.5  airflow  around  a  sphere  are  presented  in  Fig.  1.12.  If  the  heat  release 
of  a  sufficient  intensity  was  localized  in  a  small  region  before  the  blunt  body,  the  flow  structure 
was  considerably  changed  representing  a  pattern  typical  of  the  flow  structure  in  the  presence  of 
the  structural  spike.  In  this  case,  one  can  see  a  departure  of  a  normal  shock  from  the  sphere  and 
its  transformation  into  the  oblique  shock,  thus,  reducing  drag.  The  prolate  subsonic  recirculation 
zone  was  predicted  (see  Fig.l.l2b)  before  the  body.  This  phenomenon  attributed  to  the  optical 
discharge  has  been  described  in  [4]. 

Then,  the  analytical  study  using  the  approach  proposed  by  Guvernuk  [13]  was  carried  out 
allowing  to  conclude  that  the  drag  reduction  effect  is  achieved  mainly  due  to  the  formation  of 
the  wake  type  non-uniformity  in  the  incoming  air  stream.  The  initial  non-uniformity  can  defor¬ 
med  by  the  fuel  counterflow  injection.  In  the  presence  of  combustion,  this  type  of  non-uniformity 
can  be  reinforced  or  reduced  depending  on  the  combustion  wave  propagation  direction  and  the 
combustion  regime,  i.e.,  deflagration  or  detonation.  It  is  known  that  in  the  detonation  regime. 


(a)  (b) 

Figure  1.12:  Calculated  Mach  number  (15  isolines  ranged  between  0.0-  2.5)  distributions  around 
a  spherea)  without  and  b)  with  the  energy  release  zone  placed  before  the  body. 
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combustion  products  move  in  the  direction  of  the  propagating  combustion  wave.  I  n  the  deflag¬ 
ration  regime,  another  situation  takes  place:  gases  escaping  the  combustion  zone  are  moving 
directly  opposite.  Using  the  counterflow  gas  injection  of  combustible  gases,  it  is  difficult  to  ex¬ 
pect  that  the  detonation  regime  could  be  realized.  If,  in  the  deflagration  regime,  counterflow 
combustion  is  realized  by  forming  the  flame  propagating  to  the  body  surface,  the  "flame"  spike 
effect  will  rei nforce  the  effect  of  hydrogen  i nj ecti on . 

The  term  "deflagration"  summarizes  flame  propagation  regimes  corresponding  to  a  low  branch 
of  the  Hugoniot  curve.  Recently  [7],  the  enhanced  turbulent  deflagration  mode  was  numeri¬ 
cally  predicted  in  the  case  of  a  stoichiometric  hydrogen/air  mixture  under  atmospheric  condi¬ 
tions  when  ignition  was  initiated  by  a  powerful  electrical  discharge  in  a  duct  constructed  as  a 
'two-room  compartment",  where  two  volumes,  donor  and  receiver,  were  connected  by  a  narrow 
passage.  The  mathematical  model  and  modifications  of  the  Kl  VA-3  computer  code  [6]  used  for 
numerical  simulation  are  described  in  [7].  In  order  to  account  the  effect  of  turbulence/chemistry 
interactions,  thePaSR  (Partially  Stirred  Reactor)  model  in  a  generalized  form  has  been  incorpo¬ 
rated  into  a  framework  of  the  RNG  A;  -  e  turbulence  model  for  reacting  compressible  gases,  to 
simulation  hydrogen  counterflow  combustion.  The  brief  model  description  is  given  below. 

1.4.1  Model  of  turbulent  combustion 

M  ost  methods  i  n  react!  ng  flow  model  i  ng  i  nvol  ve  numeri  cal  sol  uti  ons  of  the  Reynol  ds-averaged 
Navier-  Stokes  (RANS)  equations.  More  sophisticated,  but  more  promising  approaches  make  use 
of  the  probability  density  function  (PDF)  technique  coupled  with  the  large  eddy  simulation  (LES) 
models.  Once  the  PDF  is  known,  the  mean  chemical  reaction  rates  can  be  treated  in  a  closed 
form,  irrespective  of  complexity  of  a  chemical  mechanism.  Some  methods  use  a  prescribed  form 
of  the  PDF  parameterized  by  its  lower  statistic  moments,  in  others,  the  function  is  defined  by 
solving  the  evolutional  single-point  PDF  equation.  Both  approaches  become  problematic  when 
applied  to  complex  chemistry.  With  the  assumed  PDF  method,  the  choice  of  the  PDF  form  is 
quite  arbitrary,  and,  as  the  number  of  species  grows,  this  technique  leads  to  an  unlimited  num¬ 
ber  of  low-moment  equations  which  must  be  closed  and  solved.  For  thedirect  PDF  approach,  the 
problem  arises  when  the  unclosed  micro-  mixing  terms  are  modeled. 

Application  of  another  useful  idea  to  turbulence  combustion  modeling  assumes  that  under 
high-intensity  conditions,  turbulence  exerts  the  main  impact  on  the  mechanism  of  turbulence 
combustion,  but  the  influence  of  finite-rate  chemistry  is  never  negligible.  Theformation  of  an  ar¬ 
chipelago  of  unburnt  gas  pockets  dueto  the  distorted  flame  front  reconnections  can  be  regarded 
as  a  main  consequence  of  the  model.  I  n  this  case,  thePaSR  (Partially  Stirred  Reactor)  model  has 
been  generalized  to  account  for  the  effect  of  mixture  imperfections  on  chemical  reaction  rates. 

The  model  distinguishes  (see  Fig.l.lBa)  between  the  concentration  (in  mean  molar  density) 
an  the  reactor  exit,  the  concentrations  in  the  reaction  zone,  c,  and  in  the  feed,  (9.  To  il¬ 
lustrate  the  relation  between  rates  of  the  processes  I -I  I  in  the  reactor,  let  us  analyze  the  PaSR 
differential-  algebraic  problem 


Tq  Tfnix  Tc 

where  T  is  the  time  step,  Tc  is  the  chemical  reaction  time,  and  Tmix  is  the  micro-mixing  time. 

The  equations  above  (represented  graphically  in  Fig.l.lBb),  express  the  idea  that  combustion 
in  the  reactor  is  a  sequential  process  where  mixing  is  followed  by  chemical  reactions,  and  that 
in  any  moment,  the  rates  of  the  individual  steps  are  equal.  After  algebraic  manipulations,  one 
can  yield  the  analytical  solutions  of  the  I  inear  problem  (1.1): 
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Figure  1.13:  a)  The  concept  of  a  partially  stirred  reactor,  PaSR.  The  reaction  zone  is  lined;  b) 
Process  rates  diagram  for  the  PaSR. 


The  sol  uti  on  (1.2)  shows  that  the  turbul  ent  combust!  on  ti  me  i  s  a  sum  of  the  mixing  and  reaction  times, 
if  the  process  is  expressed  in  terms  of  the  reactor  output  parameters.  Thus,  the  model  is  similar 
to  turbulent  dissipation  (or  eddy  break-up)  approach,  however,  accounting  for  finite-rate  che¬ 
mistry.  Finally,  the  sol  uti  on  (1.2)  can  be  represented  as 


—  c° 


T 


1  rtl 

Tc  ^  Tc  Tmix 


(1.3) 


where  K  =  Tc/irmix  +  Tc),  and  Hm  is  a  harmonic  mean. 

The  model  generalization  is  equivalent  to  a  substitution  of  the  equations  (1.1)  for  a  set  of  the 
general  PaSR  rate  balances 


/r-(c),  (1.4) 

'^mix 

where 

fr{c)  =  {f”  —  ul)ujr{c)  =  termf  —  termf  (1.5) 

is  the  Arrhenius  chemical  source  term  separated  into  production,  ternf,  and  destruction,  termf, 
rates,  ly'f  and  are  the  stoichiometric  coefficients  of  the  backward  and  forward  stages,  respecti¬ 
vely.  Species  indices  are  omitted  for  simplification. 

The  r-reaction  progress  variable ^^(c)  is  calculated  from  the  mass-action  law 

Ns  Ns 

Uric) = 4iT)  n  (q)"”  -  kiiT)  n 

i=l  i=l 

where  T  is  the  temperature,  fc/  and  are  the  rate  coefficients  for  the  forward  and  backward 
stages  of  the  reaction,  and  Ns  is  the  total  number  of  species  in  the  mixture. 

The  first  equation  of  the  system  (1.4)  represents  the  chemical  step  of  the  operator  splitting 
procedure  applied  to  the  species  mass  balance.  The  second  one  is  the  time  dependent  species 
mass  conservation  equation  analogous  to  that  in  the  PSR  model,  when  the  residence  time  is 
replaced  by  the  micro-mixing  time.  From  matching  the  systems  (1.1)  and  (1.4)  follows  that  they 
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are  analogous,  if  the  second  equation  in  the  system  (1.4)  is  taken  in  a  steady-state  form.  Finally, 
one  can  get  the  general  rate  expression 


c]  -Cs 


=  fr{cl)  ■  K  = 


c°  +  termfr  +  termfr'^ 
u _ /  '-sJr  \ 

Hm[  ©  ’  j.  0  /  ^ 

^  ci+termrT  terrurr’i^ 


(1.6) 


containing  no  reaction  zone  parameters  c,  which  cannot  be  resolved  on  a  computational  grid, 
but  replacing  their  effect  with  the  rate  multiplier  k  =  Tc/(tc  -t  =  l/2Hm{t,Tmix)- 

The  superscri  pt  °  denotes  the  val  ues  at  the  start  of  the  ti  me  i  ntegrati  on,  whi  I  e  the  subscri  pt  * 
refers  to  the  "reference"  speci  es  whi  ch  are  equi  val  ent  to  the  "I  i  miti  ng"  speci  es  of  the  M  agnussen- 
Hjertager  E DC-model. 

The  formal  i  ntroducti  on  to  the  concept  of  "reference"  speci  es  fad  I  itati  ng  calcul  ations  of  reac- 
tion  rates  in  the  multi-stage  chemical  mechanism  can  be  found  in  [14].  If  Tmix  0,  the  model 
reduces  to  the  quasi-laminar  approach  with  =  c  in  the  reaction  rate  terms.  Similarly  to  Eq. 
(1.3),  the  arguments  of  the  harmonic  mean  function  in  Eq.  (1.6)  are  the  reaction  rate  and  the 
eddy  dissipation  rate,  respectively.  In  the  limit  of  fast  chemistry,  i.e.,  ternf  ■  c°  =  termf  ■  c],  the 
dissipation  ratein  Eq.(1.6)  reproduces  the  classic  eddy  break-up  rate,  ~  (4  -cfi/Tmix- 


1.4.2  Definition  of  micro-mixing  time 

IftheRNG  k-e  model  is  employed,  the  turbulent  viscosity  related  to  k,  the  turbulent  kinetic 
energy,  ande,  the  dissipation  rate  of  k,  is  given  by  the  general  expression. 

I^t  =  =  p{vi  +  (1.7) 

V 

where  is  the  "standard"  k-e  value  of  the  kinematic  viscosity. 

The  above  expression  represents  contributions  of  three  mechanisms  of  diffusion  transport, 
molecular,  average  turbulent,  and  the  additional  one  responsible  for  micro-mixing.  The  model  is 
assumed  to  be  valid  across  a  full  range  of  flow  conditions  from  low  to  high  Reynolds  numbers,  if 
k  and  e  are  deter  mi  ned  from  the  general  i  zed  transport  equati  ons, 

P-pe  +  V  ■[{^)Vk]  (1.8) 

c,,,{pe)V  •  U  +  -  c,,pe)  +  V  •  [(;^)Ve], 

where  P  is  the  production  term,  given  by 

P  =  Pt[S  -\{V  ■\jf]-‘^pkV  -  v, 
and  S  is  the  mean  strain  rate,  defined  as 

1  ^  ^  2 
2^  dxi  dxj  ’ 

Provided  that  thefirst  two  terms  in  (1.7)  contributed  to  the  conventional  diffusion  transport, 
the  geometrical  mean  term  in  the  expression  can  be  used  to  define  the  characteristic  time  scale 
for  micro-mixing  in  the  turbulence/chemistry  interaction  model,  if  written  in  a  form, 

'^mix  ~  jk  =  2Cn  ^  T/;  ~  T/;,  (1-9) 
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where  =0.09,  Tk  =  {vileYl‘^  is  the  Kolmogorov  time. 

Thus,  the  sequence  of  diffusion  processes  is  represented  by  different,  gradient  and  non¬ 
gradient,  approaches.  Eg.  (1.9)  can  be  rewritten  using  the  eddy  break-up  time  definition 

'^mix  ^  Tk  ~  {C/j  lRetf‘^kle  (1.10) 

This  leads  to  a  correct  dependence  of  Tmix  on  the  Ret=k't" /k'l,  but  its  actual  value  has  to  be  opti¬ 
mized.  For  a  typical  ~  10^  case,  the  value  is  about  0.01A:/e  used  in  the  modeling.  If  smallest 
eddies  are  enlarged  to  be  resolved  by  a  computer  grid,  i.e., 

poo  poo 

has  =  /  dkk-^E{k)/  dkE{k)^A 

Jl/A  Jl/A 

POO 

usgs  =  /  dkk-^[k-^E{k)Y/^ 

Jl/A 

another  definition  of  the  micro-mixing  time  has  been  used. 

Tmix  ~  A2/3e-V3  (1.11) 

where  A  is  the  minimal  scale  resolved  on  the  grid,  E{k)  =  CS‘l'^k~^l'^  is  the  Kolmogorov  energy 
spectrum.  Theformula  (1.11)  renders  the  model  a  quality  of  the  SGS  (sub-grid  scale)  approach. 

1.4.3  Results  and  discussion 

To  account  for  the  hydrogen/oxygen  distribution  and  combustion  after  hydrogen  counterflow 
injection,  both  mixing  and  ignition/combustion  processes  must  be  simulated.  To  simulate  the 
effect  of  combustion,  the  reaction  mechanism  of  hydrogen  oxidation  in  air  has  been  constructed. 
M  ost  of  the  react!  ons  have  been  studi  ed  experi  mental  I  y,  and  the  mechani  sm  has  previ  ousi  y  been 
used  for  hydrogen  propulsion  modeling  [9].  The  reaction-rate  parameters  in  the  form  of  kj  = 
exp{-Eaj I RqT)  are  taken  from  different  sources  and  the  electronic  copy  of  the  mechanism 
can  be  found  on  our  web  site  [10]. 

Using  the  reaction  mechanism  developed,  the  calculated  ignition  delay  times  for  stoichio¬ 
metric  hydrogen  air  mixtures  at  different  initial  pressures  and  temperatures  were  calculated 
and  compared  with  the  shock-tube  experiments.  Good  agreement  has  been  achieved  as  illustra¬ 
ted  in  Fig.  1.14.  The  main  goal  in  using  the  detailed  mechanism  of  hydrogen  combustion  is  to 


Figure  1.14:  Validation  of  H2/air  ignition  chemistry.  Ignition  promotion  and  moderation  of  its 
dependence  on  pressure  for  the  stoichiometric  mixture  by  adding  a  small  amount  of  silane. 


13 


correctly  reproduce  rather  ill-ignition  conditions  above  the  so-called  extended  second  explosion 
limitfor  hydrogen/air  mixtures.  Thiseffect  iswell  pronounced  in  the  data  presented  in  Fig.  1.14. 
Duetothis  effect,  the  ignition  delay  for  higher  pressure  which  is  shorter  than  the  ignition  delay 
for  lower  pressure  at  temperatures  higher  than  the  explosion  limit  temperature,  becomes  longer, 
if  the  initial  temperatures  are  lower  than  the  limit.  This  is  why,  one  can  seethe  intersection  of 
two  experimental  (for  different  pressures)  curves  presented  in  Fig.  1.14.  This  phenomenon  is 
likely  responsi bl e for  deterioration  of  drag  reduction  efficiency  observed  in  the  experiments  (see 
Fig.  1.7),  if  larger  amounts  of  hydrogen  were  injected.  As  a  result,  this  can  lead  to  the  incre¬ 
ased  pressure  level  in  the  ignition  zone,  thus,  worsening  auto-ignition  conditions.  To  moderate 
the  ignition  delay  dependence  on  pressure,  it  has  been  assumed  in  [12]  to  make  use  of  si  lane  as 
the  ignition  promoter.  The  measure  facilitates  the  ignition  process  which,  in  this  case,  becomes 
nearly  hypergolic. 

The  results  of  2-D  counterflow  combustion  simulations  using  the  Kl  VA3  computer  code  mo¬ 
dified  to  account  for  detailed  chemistry  hydrogen/air  combustion  are  presented  in  Fig.  1.15.  I  n 


■ 


(a)  t=0.4  ms 


(b)  t=0.5  ms 


(c)  t=0.7  ms 


(d)  t=1.6  ms 


(e)  t=2.5  ms 


(f)  t=0.8  ms 


(g)  t=1.0  ms 


(h)  t=1.0  ms 


(i)t=1.0ms 


(j)t=1.0  ms 


Figure  1.15:  Time  development  of  couterflow  H2/air  combustion  for  the  pipe-spike.  Plots  a)-e) 
correspond  to  pressure  distributions  at  different  moments;  plots  f)-g)-  to  temperature  distribu¬ 
tions;  plots  h),  i)  and])-  to  distributions  of  hydrogen  peroxide,  water  vapor  and  atomic  hydrogen, 
respectively. 
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the  axi symmetric  simulations,  a  part  of  a  cylinder  body  of  a  size  2x2  cm  with  a  structural  spike 
of  a  length  2  cm  was  placed  in  a  computation  region  of  a  size  8x20  cm.  The  structural  spike  was 
a  pipe  of  the  inner  radius  0.25  cm,  and  the  outer  radius  0.5  cm.  Hydrogen  has  been  injected 
through  the  pipe  with  a  speed  of  50  m/s.  The  free  air  stream  velocity  corresponded  to  M=2.0 
and  Po=l-o  atm.  Hydrogen  ignition  has  been  initiated  by  two  hot  spots:  the  first  one  was  pla¬ 
ced  inside  the  pipe  (for  a  short  time),  and  the  second  one  -  0.5  cm  below  the  spike  end  in  the 
recirculation  zone.  Transient  calculations  continued  till  the  steady-state  was  achieved  (but  not 
longer  than  10  ms).  All  this  time,  hydrogen  injection  continued,  but  the  external  hot  spot  was 
supported  by  the  energy  source  in  the  course  of  1  ms.  Within  the  ignition  window,  the  specific 
internal  energy  in  the  specified  mesh  cells  was  increased  on  each  timestep.  If  the  temperature 
in  the  ignition  cell  was  i n  excess  of  2000  K  before  the  end  of  the  ignition  window,  then,  in  this 
point,  the  energy  deposition  was  terminated. 

The  block-structured  mesh  KIVA-3  code  used,  in  the  modeling,  might  be  considered  of  only  a 
moderate  accuracy  for  aerodynamic  calculations,  but  such  an  accuracy  seems  to  be  sufficient  to 
si  mulatethe  effects  of  counterflow  combustion  on  the  supersonic  flow  structure.  From  the  results 
presented  in  Fig.  1.15,  one  can  conclude  that  the  typical  shock  structure  (see  Fig.  1.15a)  formed 
in  front  of  the  blunt  body  at  the  initial  moment  was  smoothed  later  on  while  combust!  on  process 
was  developing  (see  plots  Fig.  1.15b-c).  This  prediction  is  in  agreement  with  the  experimental 
observation.  From  the  temperature  distributions  (see  plots  in  Fig.  1.15f-g),  one  can  see  that  the 
flame  was  stabilized  in  the  frontal  recirculation  zone  formed  in  the  presence  of  the  structural 


(a)  t=0.2  ms  (b)  t=0.6  ms  (c)  t=0.7  ms  (d)  t=0.8  ms  (e)  t=1.0  ms 


(f)t=1.2ms  (g)t=1.4ms  (h)t=1.7ms  (i)t=2.0ms  (j)t=2.5ms 


Figure  1.16:  Time  development  of  couterflow  H2/air  combustion  for  the  conical  tip  spike.  Plots 
a)-i )  correspond  to  temperature  distri  buti  ons  at  different  moments;  the  plot  j )-  to  the  distri  bution 
of  water  vapor. 
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spike.  The  distributions  of  intermediate  (hydrogen  peroxide)  species  and  main  (water  vapor) 
combustion  products  as  well  as  a  high  concentration  of  H-  radicals  show  that  the  heat  release 
takes  place  only  in  this  region,  effectively  forming  the  flame  spike  structure  of  the  counterflow 
combustion.  The  temperature  in  the  reaction  zone  was  predicted  to  be  only  ~  1300  K.  The 
predicted  flame  shape  and  position  testifies  to  the  fact  that  combustion  takes  place  in  the  lean 
mixture,  and  the  deflagration  regime  is  realized  with  the  flame  front  moving  towards  to  the  body 
surface.  In  this  case,  the  reaction  products  are  moving  in  the  same  direction  as  the  injected 
gases,  thus,  increasing  the  positive  effect  of  injection.  As  the  high  level  of  temperature  in  the 
reaction  zone  was  not  achieved  at  modeling  conditions,  special  measures  should  be  taken  to 
enhance  and  stabilize  combust!  on.  To  stabilize  the  flame  in  the  experi  ments,  a  conical  tip  on  the 
spike  end  was  used.  Thus,  the  spike  geometry  was  modified  in  the  modeling  to  account  for  the 
presence  of  the  tip,  and  the  computation  results  are  presented  in  Fig.  1.16. 

In  this  modeling,  the  spike  length  was  increased  by  1.5  cm,  from  which  0.5  cm  was  a  shear 
of  the  conical  tip.  The  pipe  inner  diameter  was  reduced  to  0.2  cm,  that  allowed  to  specify  a 
higher  injection  velocity  (200  m/s)  without  substantial  increasing  of  the  hydrogen  mass  flow 
rate.  The  simulation  scenario  was  also  changed:  hydrogen  injection  and  ignition  (see  Fig.  1.16a) 
took  place  in  the  quiescent  atmosphere.  Then,  at  the  moment  t=0.5  ms  the  air  was  set  in  motion 
corresponding  to  M  =2.0  and  the  total  pressure  Po=  1.0  atm.  The  air  stream  parameters  were 
calculated  using  the  isentropic  flow  relations  for  a  perfect  gas.  The  air  stream  formed  the  region 
of  elevated  pressure  and  temperature  which  interacted  with  the  developed  hydrogen  flame  as 
can  be  seen  in  Fig.l.l6b-c.  As  the  result  of  the  interact!  on,  the  flame  front  was  transform  into  a 
typical  form  of  a  "butterfly",  the  shape  of  which  was  observed  in  the  experi  ments.  Then,  incoming 
air  flow  blew  the  flame  almost  away,  and  combustion  continued  only  in  the  pipe.  Later  on  (see 
Fig.  1.16d-e),  the  hydrogen  combustion  in  the  supersonic  stream  was  re-established,  and  the 
reaction  zone  structure  rebuilt  the  configuration  typical  of  the  developed  flame  spike  (see  Fig. 
1.16f-i).  Water  vapor  distribution  in  Fig.  1.16j  illustrates  an  important  role  of  the  tip  in  the 
flame  stabilization.  The  higher  temperature,  1420  K,  was  predicted  in  the  react! on  zone. 

The  simulation  results,  when  the  hydrogen  injection  velocity  was  increased  iup  to  250  m/s 
are  presented  in  Fig.  1.17.  As  a  result,  the  combustion  zone  became  more  extended  (see  Fig. 
1.17a),  but  a  character  of  the  flame  front  interaction  with  the  air  stream  remain  unchanged,  e.g., 
at  first,  the  flame  "butterfly"  structure  was  formed,  then,  the  flame  was  almost  extinguished, 
and  regenerated  once  again  (see  Fig.  1.17b-e).  The  temperature  maximum  in  the  reaction  zone 
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(a)  t=0.6  ms  (b)  t=0.7  ms  (c)  t=0.8  ms  (d)  t=1.0  ms  (e)  t=2.5  ms 

Figure  1.17:  Time  development  of  couterflow  H2/air  combustion  for  the  conical  tip  spike.  Hydro¬ 
gen  njection  velocity  is  equal  to  250  m/s.  Plots  a)-i)  correspond  to  temperature  distributions  at 
different  moments. 
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reached  1840  K,  and  combustion  was  well  stabilized  in  the  recirculation  zone  behind  the  tip.  One 
can  expect  than  this  tendency  could  be  extrapol ated  to  the  higher  hydrogen  injection  velocities. 
However,  it  requires  further  reduction  of  the  pipe  inner  size  that  arises  a  problem  of  the  mesh 
resolution  in  simulations  of  a  moderate  accuracy. 

As  a  concei  vabi  e  measure  of  temperature  i  ncrease  i  n  the  combust!  on  zone  coul  d  be  consi  dered 
an  addition  of  boron  to  hydrogen  injected  through  the  spike.  The  boron  effect  on  the  hydrogen 
combustion  has  been  studied  using  the  time-dependent  reaction  code  Senkin  of  the  Chemkin 
library,  and  the  results  of  computations  are  presented  in  Fig.  1.18.  The  extended  chemical  me- 


Auto-ignition  of  Hj-Boron-air  Mixture 


(a) 


Auto-ignition  of  Hj-Boron-air  Mixture 


(b) 


Auto-ignition  of  Hj-Boron-air  Mixture 

(P„=1.0  bar,  T„=800  K,  H2:B2=3:1) 


(C) 


Figure  1.18:  a)  Selected  concentrations  of  combustion  products  vs  time  histories,  b)  selected  con¬ 
centrations  of  combustion  products  vs  time  hi  stories,  c)  illustration  of  reduced  ignition  delay  and 
increased  heat  release  when  compared  with  parameters  of  similar  stoichiometric  hydrogen/air 
mixture. 


chanism  accounting  for  the  presence  of  boron  vapor  in  the  mixture  can  be  found  on  our  web  site 
[10].  The  comprehensive  mechanism  of  hydrogen/boron  combustion  consists  of  64  elementary 
reactions  between  22  species  including  B,  B2,  BO,  BO2,  B2O,  B2O2,  B2O3,  B203(L),  B(L),  L  deno¬ 
tes  the  liquid  phase.  The  auto-  ignition  of  H2/boron  mixture  (3:1,  in  mole  numbers)  was  studied 
with!  n  the  i  niti  al  temperature  range  of  700-  900  K ,  and  the  pressure  range  of  1-  5  bar.  I  n  the  F  i  g. 
1.18a-b,  one  can  see  that  the  concentration  of  B2O3  in  the  gas  phase  exceeds  that  of  0203(1 )  at 
the  equilibrium  state.  However,  there  is  a  moment  in  the  development  of  the  auto-ignition  pro¬ 
cess  development,  when  these  phases  are  predicted  in  comparable  quantities.  Such  a  presence 
of  the  liquid  phase  in  combustion  products  can  reinforce  the  effect  of  the  "flame"  spike. 

The  results  presented  in  Fig.  1.18c  illustrate  that  combustion  of  H2/boron  composition  has 
essential  advantages,  if  compared  with  combustion  of  H2/air  mixtures:  theignition  delays  are  by 
a  order  of  magnitude  shorter,  and  the  temperature  of  combustion  products  is  by  1000  K  higher. 
It  is  in  order  to  remark  that  the  use  of  silane  daes  not  affect  the  heat  release  in  the  reactions. 
In  contrast  to  H2/air  mixtures,  the  increase  in  the  initial  pressure  up  to  5  bar  does  not  affect 
much  the  auto-ignition  process.  All  these  factors  are  beneficial  of  the  application  of  theH2/boron 
compositions  to  increase  the  positive  effect  of  the  counterflow  combustion  on  drag  reduction  of 
the  blunt  bodies  in  supersonic  streams. 

The  theoreti  cal  and  computati  onal  anal  ysi  s  has  been  car  ri  ed  out  at  the  department  of  Ther  mo- 
and  Fluid  Dynamicsof  the  Chalmers  University  of  Technology. 
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1.5  Conclusions 


•  The  theoretical  and  experimental  study  of  counterflow  combustion  in  supersonic  streams 
confirms  the  fact  that  the  "flame"  spike  can  be  considered  as  an  effective  means  to  reduce 
the  pressure  drag  of  blunt  bodies.  Direct  measurements  of  pressure  drag  have  been  carried 
out  in  the  supersonic,  M  =2. 0-2.5  air  streams  when  the  "flame"  spike  has  been  formed  with 
the  help  of  counterflow  combustion  of  hydrogen  injected  through  thesmall  pipe  playing  role 
of  the  structural  spike.  Thespikewith  a  conical  tip  having  better  flame  stabilization  capa¬ 
bilities  has  been  also  studied.  The  main  mechanism  of  drag  reduction  has  been  analyzed 
with  the  help  of  numerical  simulations. 

In  experiments,  the  regions  of  stable  counterflow  combustion  of  hydrogen  in  the  d  (the 
blunt  body  diameter)  -  L  (the  dimensionless  length  of  the  spike)  andG  (the  dimensionless 
hydrogen  mass  flow  rate)  -  Po  (free  stream  stagnation  pressure)  planes  have  been  deter¬ 
mined.  The  main  elements  of  the  flow  structure  around  the  blunt  body  in  the  presence  of 
inert  and  reactive  injection  have  been  resolved  and  measured. 

•  With  the  help  of  theoretical  analysis,  it  has  been  established  that  the  mechanism  of  drag 
reduction  is  mainly  attributed  to  formation  of  the  wake-type  non-uniformity  in  the  free 
stream  initiated  by  injected  hydrogen  and  reinforced  by  motion  of  combustion  product  be¬ 
hind  the  flame  front  propagating  downstream  in  the  deflagration  regime  of  combustion. 

When  pressure  in  the  reaction  zone  increases  (it  can  happen  with  an  increase  of  the  fuel 
mass  flow  rate  or  with  an  increase  of  free  stream  pressure),  the  hydrogen  combustion  be¬ 
comes  slower  and  less  intensive  that  results  in  reduction  of  the  positive  effect  of  the  "flame" 
spike.  To  prevent  this  consequence,  the  moderation  of  pressure  dependence  of  hydrogen 
combustion  above  the  "second  extended  limit  of  explosion"  is  proposed  using  silane  or  bo¬ 
ron.  I  n  the  case  of  boron,  one  can  expect  more  rapid  combustion  as  well  as  much  stronger 
heat  release  effect.  Both  factors  are  beneficial  of  the  drag  reduction  effect. 

•  The  results  of  this  study  illustrating  that  the  comparable  efficiency  in  drag  reduction  can 
be  achieved  in  the  same  wind  tunnel  under  the  similar  flow  conditions  by  using  two  dif¬ 
ferent  energy  deposition  methods  -  the  optical  (laser)  discharge  and  counterflow  hydrogen 
combustion.  Thus,  it  may  be  argued  that  the  phenomenon  of  drag  reduction  of  blunt  bodies 
in  supersonic  flows  at  atmospheric  conditions  is  not  plasma  specific. 

The  results  illustrate  the  efficiency  of  counterflow  combustion  in  drag  reduction  of  blunt 
bodies  in  supersonic,  M=2.0-2.5,  flows.  If  realized,  the  technique  could  help  in  removing 
the  controversial  DU  -  projectiles  from  active  applications. 

Software  developed  or  modified  as  a  part  of  this  project  can  be  delivered  if  requested. 
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